:arXiv:1103.2713vl [astro-ph.SR] 14 Mar 2011: 



Table 3. System parameters derived using the probability density functions shown in Fig. 4. R r is the volume radius of the secondary's Roche lobe (Eggleton 1983). The errors quoted 
are statistical errors of the fitting process and do not include systematic effects such as those discussed in section 3.3. 





CTCV J1300-3052 


CTCV J2354-4700 


SDSS J1152+4049 


OU Vir 


DV UMa 


XZ Eri 


SDSS 1702 


Q 


0.240 ±0.021 


0.1097 ±0.0008 


0.155 ±0.006 


0.1641 ±0.0013 


0.1778 ±0.0022 


0.118 ±0.003 


0.248 ± 0.005 


M W /M Q 


0.736 ± 0.014 


0.935 ±0.031 


0.560 ±0.028 


0.703 ±0.012 


1.098 ±0.024 


0.769 ±0.017 


0.91 ±0.03 


Rw 1 Rq 


0.01111 ± 0.00018 


0.0089 ± 0.0003 


0.0135 ± 0.0004 


0.01191 ±0.00017 


0.00703 ± 0.00028 


0.01081 ± 0.00022 


0.0092 ± 0.0004 


log g 


o.Zl it U.Uz 


O fl ±(1 f|/l 

o.Ol it U.U4 


( .yo ± u.uo 


o. lo it U.Uz 


o./o ± U.U4 


q ok -U n HQ 
o.ZD it U.uo 


o.4 ( it U.UO 


T W (K) 


11100 ±800 


14800 ± 700 


12400 ± 1400 


22300 ± 2100 


15500 ± 2400 


15300 ± 1900 


15200 ± 1200 


M t /Mq 


0.177 ±0.021 


0.101 ±0.003 


0.087 ± 0.006 


0.1157 ±0.0022 


0.196 ±0.005 


0.091 ±0.004 


0.223 ±0.010 


Rr/ Rq 


0.215 ± 0.008 


0.1463 ±0.0016 


0.142 ±0.003 


0.1634 ±0.0010 


0.2176 ±0.0018 


0.1350 ±0.0018 


0.252 ± 0.004 


0,/Rq 


0.813 ± 0.011 


0.692 ±0.008 


0.606 ± 0.010 


0.686 ± 0.004 


0.892 ± 0.006 


0.621 ±0.005 


0.945 ±0.012 


K w (kms -1 ) 


90 ±8 


51.9 ±0.6 


60 ±3 


66.4 ±0.6 


78.9 ±1.0 


53.6 ± 1.5 


94.0 ±2.2 


K r (kms" 1 ) 


372.2 ±2.5 


482 ± 6 


387 ±6 


403.0 ±2.3 


443 ±3 


452 ±3 


380 ±4 


i° 


86.3 ± 1.1 


89.26 ±0.28 


82.38 ±0.23 


79.60 ± 0.04 


82.93 ±0.10 


80.02 ±0.12 


82.55 ± 0.17 


d (pc) 


375 ± 13 


674 ± 19 


543 ± 21 


570 ± 70 


504 ± 30 


371 ± 19 


270 ± 16 


P orb (mins) 


128.0746325(14)* 


94.3923889(14)* 


97.518753(4)* 


104.696803(7) 1 


123.6278190(20) 2 


88.069667(7) 2 


144.11821(13) 3 




SDSS 1035 


SDSS 1507 


SDSS 0903 


SDSS 1227 


SDSS 1433 


SDSS 1502 


SDSS 1501 


q 


0.0571 ± 0.0010 


0.0647 ± 0.0018 


0.113 ±0.004 


0.1115 ±0.0016 


0.0661 ±0.0007 


0.1099 ±0.0007 


0.101 ±0.010 


M W /M Q 


0.835 ± 0.009 


0.892 ±0.008 


0.872 ±0.011 


0.796 ±0.018 


0.865 ± 0.005 


0.709 ±0.004 


0.767 ±0.027 


Rw/Rq 


0.00991 ± 0.00010 


0.00956 ± 0.00013 


0.00947 ± 0.00019 


0.01052 ± 0.00022 


0.00962 ± 0.00006 


0.01145 ±0.00005 


0.0107 ± 0.0003 


log g 


8.37 ±0.01 


8.45 ±0.01 


8.42 ±0.02 


8.29 ±0.02 


8.41 ±0.01 


8.17 ±0.01 


8.26 ±0.04 


T W (K) 


10000 ± 1100 


11300 ± 1000 


13300 ± 1700 


15900 ± 1400 


12700 ± 1500 


11800 ± 1200 


10800 ± 1500 


M r /M Q 


0.0475 ± 0.0012 


0.0575 ± 0.0020 


0.099 ± 0.004 


0.0889 ± 0.0025 


0.0571 ± 0.0007 


0.0781 ± 0.0008 


0.077 ±0.010 


Rr/ Rq 


0.1047 ±0.0008 


0.0969 ±0.0011 


0.1358 ±0.0020 


0.1365 ±0.0013 


0.1074 ±0.0004 


0.1241 ±0.0003 


0.122 ±0.005 


a/R Q 


0.5977 ± 0.0022 


0.5329 ± 0.0019 


0.632 ± 0.003 


0.640 ±0.005 


0.5869 ±0.0012 


0.5844 ± 0.0013 


0.588 ± 0.008 


K w (kms -1 ) 


28.5 ± 0.6 


35.1 ± 1.0 


54.6 ±2.0 


51.3 ± 0.8 


33.8 ±0.3 


50.4 ±0.4 


48 ± 5 


K r (kms -1 ) 


499.3 ± 1.5 


543.7 ± 1.2 


481.7 ± 1.9 


460 ±3 


511.1 ±0.9 


456.5 ±0.8 


470.5 ± 3.2 


i° 


83.98 ± 0.08 


83.47 ±0.12 


82.09 ±0.19 


84.29 ± 0.10 


84.36 ±0.05 


88.35 ±0.17 


82.8 ± 0.5 


d (pc) 


174 ± 12 


168 ± 12 


299 ± 14 


400 ± 13 


226 ± 12 


175 ± 11 


306 ± 21 


P orb (mins) 


82.08965(29) 4 > 6 


66.61192(6) 5 - 6 


85.065902(13) 6 


90.661019(10) 6 


78.106657(3) 6 


84.82984(7) 6 


81.85141771(28)* 



*This paper; iFeline et al. 2004a; 2 Fclinc ct al. 2004b; 3 Littlefair ct al. 2006a; 4 Littlcfair et al. 2006b; 5 Littlefair et al. 2007; 6 Littlefair et al. 2008; 
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ABSTRACT 

We present high-speed, three-colour photometry of the eclipsing cataclysmic vari- 
ables CTCV J1300-3052, CTCV J2354-4700 and SDSS J115207.00+404947.8. These 
systems have orbital periods of 128.07, 94.39 and 97.52 minutes respectively, placing 
all three systems below the observed "period gap" for cataclysmic variables. For each 
system we determine the system parameters by fitting a parameterised model to the 
observed eclipse light curve by x 2 minimisation. 

We also present an updated analysis of all other eclipsing systems previously 
analysed by our group. The updated analysis utilises Markov Chain Monte Carlo 
techniques which enable us to arrive confidently at the best fits for each system with 
more robust determinations of our errors. A new bright spot model is also adopted, that 
allows better modelling of bright-spot dominated systems. In addition, we correct a 
bug in the old code which resulted in the white dwarf radius being underestimated, and 
consequently both the white dwarf and donor mass being overestimated. New donor 
masses are generally between 1 and 2er of those originally published, with the exception 
of SDSS 1502 (-2.9cr, AM r = -0.012M Q ) and DVUMa (+6.1<r, AM r = +O.O39M ). 
We note that the donor mass of SDSS 1501 has been revised upwards by O.O24M 
(+1.9er). This system was previously identified as having evolved passed the minimum 
orbital period for cataclysmic variables, but the new mass determination suggests 
otherwise. Our new analysis confirms that SDSS 1035 and SDSS 1433 have evolved 
past the period minimum for cataclysmic variables, corroborating our earlier studies. 

We find that the radii of donor stars are oversized when compared to theoretical 
models, by approximately 10 percent. We show that this can be explained by invoking 
either enhanced angular momentum loss, or by taking into account the effects of star 
spots. We are unable to favour one cause over the other, as we lack enough precise 
mass determinations for systems with orbital periods between 100 and 130 minutes, 
where evolutionary tracks begin to diverge significantly. 

We also find a strong tendency towards high white dwarf masses within our sample, 
and no evidence for any He-core white dwarfs. The dominance of high mass white 
dwarfs implies that erosion of the white dwarf during the nova outburst must be 
negligible, or that not all of the mass accreted is ejected during nova cycles, resulting 
in the white dwarf growing in mass. 

Key words: binaries: close - binaries: eclipsing - stars: dwarf novae - stars: low mass, 
brown dwarfs - stars: novae, cataclysmic variables - stars: evolution. 



1 INTRODUCTION 

Cataclysmic variable stars (CVs) are a class of interacting 
binary system undergoing mass transfer from a Roche-lobe 
filling secondary to a white dwarf primary, usually via a gas 
* E-mail:chris. savoury@sheffield.ac.uk stream and accretion disc. A bright spot is formed where 
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the gas stream collides with the edge of the accretion disc, 
often resulting in an 'orbital hump' in the light curve at 
phases 0.6 — 1.0 due to the area of enhanced emission ro- 
tating i nto our l i ne of sigh t. For an exc ellent overview of 
CVs, see I Warned 1 19951 ) and lHellierl (|200lh . The light curves 
of eclipsing CVs can be quite complex, with the accretion 
disc, white dwarf and bright spot all being eclipsed in rapid 
succession. When observed with time resolutions of the or- 
der of a few seconds, this eclipse structure allows the sys- 
tem parameters to be determined to a high degree of pre - 
cision with relatively few assumptions (|Wood et al.lll98a) . 
Over the last eight years our grou p has used the high -speed, 
three-colour camera ultracam |Dhillon et al.ll2007t ) to ob- 
tain such time-resolution. The ability to image in three dif- 
ferent wave-bands simultaneously makes ultracam an ideal 
tool to study the complex, highly variable light curves of 
CVs. Using ultracam data, we have obtained system pa- 
rameters for several short period systems (e.g. Feline et al. 
2004a, 2004b; Littlefair et al. 2006a, 2007, 2008), including 
the first system accreting from a sub-stellar donor (Littlefair 
et al. 2006b). 

Despite extensive study over recent decades, there are 
still several outstanding issues with evolutionary theories of 
CVs that have wide ranging implications for all close binary 
systems. The secular evolution of CVs is driven by angu- 
lar momentum losses from the binary orbit. In the standard 
model, systems with orbital periods below ~130 minutes are 
thought to lose angular momentum via gravitational radi- 
ation. Angular momentum losses sustain mass transfer and 
subsequently drive the system to shorter orbital periods, un- 
til the point wh ere the donor star becomes degenerate (e.g. 
|Paczvnski|[l98ll) . Here, the donor star is driven out of ther- 
mal equilibrium and begins to expand in response to mass- 
loss, driving the system to longer orbital periods. We there- 
fore expect to observe a period cut-off around P or t — 65 — 70 
minutes, dubbed the "period minimum", in addition to a 
build up of systems at this minimum period (the "period 
spike" ; Kolb & Bar affe 1999). The p eriod spike has recently 
been identified by iGansicke et al.1 (|2009l ). whose study of 
SDSS CVs found an accumulation of systems with orbital 
periods between 80 and 86 minutes. This is significantly 
longer than expected. A larger than expected orbital period 
implies that the orbital separation is larger than expected, 
nd thus the radii of the donor star must also be larger than 
expected in order to remain Roch e-lobe filling. Recent ob- 
servations bv lLittlefair et all |2008t ) support this, suggesting 
that the donor stars in short period CVs are roughly 10 per- 
cent l arger than predicted by the models of iKolb &: Baraffe] 
(1999). The reason why the donor stars appear oversized 
remains uncertain. Possible explanations in clude some form 
of enhanced angular momentum loss (e.g . iPattersonl 1 19981 ; 
IKolb fc BaraHl 19991 ; rWillems et al.1 [2005! ) which would in- 
crease mass loss and drive the donor stars further from ther- 
mal equilibrium, or missing stellar physics in the form of 
magn etic activity coupled with the effects of rapid rotation 
(e.g. IChabrier et ai]|2007l ). 

One way to determine why the donor stars appear 
oversized is to compare the shape of the observed donor 
mass-period relationship (A/2 — Port), and by implica- 
tion, mass-radius (M2 — R2) relationship, to the models of 
IKolb fc Baraffel (| 1999T ) . calculated with enhanced angular 
momentum loss or modified stellar physics. These models, 



in principle, make different predictions for the shape of the 
mass-period relationship and the position of the period min- 
imum. Both the shape of the mass-period relationship and 
the position of the period minimum are dependent on the 
ratio « = tm/tkh, where tm and tkh are the mass- loss and 
thermal timescales of the donor star respectively. Initially, 
kS>1, and the donor is able to contract in response to mass- 
loss. As the system evolves to shorter orbital periods both 
timescales increase, although the thermal timescale increases 
much faster than the mass-loss timescale. This results in k 
decreasing with orbital period. When the two timescales be- 
come comparable, the donor is unable to contract rapidly 
enough to maintain thermal equilibrium and becomes over- 
sized for a given mass. Since donor expansion does not occur 
until the thermal and mass-loss timescales become compa- 
rable, if enhanced angular momentum loss is responsible for 
the oversized CV donors, the systems immediately below the 
period gap would not be expected to be far from thermal 
equilibrium. In contrast, star spots would inhibit the con- 
vective processes in all CV donors below the period gap 
(assuming of course that spot properties are similar at all 
masses). Models that include the effects of enhanced angu- 
lar momentum loss and star spot coverage therefore begin to 
diverge significantly in M and AI2 at orbital periods of 100 
minutes. We can distinguish between these models if we have 
a sample of CVs that covers a wide range of orbital periods 
and whose component masses and radii are known to a high 
degree of precision (e.g. <tA/2~ O.OO5M ). Unfortunately we 
lack enough precise mass-radii determinations for systems 
with orbital periods between 95 and 130 minutes. To over- 
come this shortage we observed eclipses of three CVs below 
the period gap: CTCV J1300-3052, CTCV J2354-4700 and 
SDSS Jl 15207.00+404947.8 (hereafter CTCV 1300, CTCV 
2354 and SDSS 1152). 

CTCV 2354 and CTCV 1300 we re discovered as part 
of the Calan-Tololo Survey follow up (|Tappert et alj|2004l ). 
During the follow up, both systems were found to be eclips- 
ing with orbital periods of 94.4 and 128.1 minutes, respec- 
tively. Basic, non-time resolved, spectroscopic data was ob- 
tained for each system. The spectrum of CTCV 1300 showed 
features typical of the three main components in CVs: strong 
emission lines from the accretion disc, broad, shallow absorp- 
tion features from the white dwarf and red continuum and 
absorption bands from the donor. CTCV 2354 was found to 
contain strong emission lines of H and He, generally typical 
of a dwarf nova in quiescence. 

SDSS 1152 was identified as a CV by ISzkodv et all 
The system shows broad, double-peaked, Balmer 
emission lines, which are characteristic of a high-inclination 
accret ing binary. Follow up work by ISouthworth et all 
|2010l) found the system to have an orbital period of 97.5 
minutes. 

In this paper we present ultracam light curves of 
CTCV 1300 {u'g'r'i'), CTCV 2354 (u'g'r') and SDSS 1152 
(u'g'r'), and in each case attempt to determine the sys- 
tem parameters via light curve modelling. In addition, we 
also present an updated analysis of all eclipsing systems 
previously published by our group: OU Vir (Feline et al. 
2004a), XZ Eri and DV UMa (Feline et al. 2004b), SDSS 
J1702+3229 (Littlefair et al. 2006a), SDSS J1035+0551 
fLittlefair et a l 200 6 b, 20 08), SDSS J150722+523039 
jLittlefair et al.1 l2007l . |200S| ). SDSS J0903+3300, SDSS 
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J1227+5139, SPSS J1433+1011, SPSS J1501+5501 and 
SPSS J1502+3334 (|Littlefair et alJl200Sh . Our primary rea- 
son for doing so was the introduction of a new analysis util- 
ising Markov Chain Monte Carlo (MCMC) techniques and 
an updated bright spot model. The MCMC analysis is more 
reliable at converging to a best fit than the downhill simplex 
algorithm used previously, while the new bright spot model 
allows for a more realistic modelling of bright spot domi- 
nated systems (e.g. CTCV 1300, PV UMa, SPSS 1702) and 
should thus provide more accurate values of the mass ratio, 
q. While implementing these changes, we also discovered a 
bug in the code previously used to bin the light curves (see 
e.g. Littlefair et al. 2006a for details of the original code) 
which resulted in the white dwarf radius being underesti- 
mated, and consequently, the white dwarf and donor mass 
being overestimated. Full details are provided in section [3731 



Table 2. Orbital cphemcrides. 



2 OBSERVATIONS 

In Table [1] we present details of the observations used to 
analyse CTCV 1300, CTCV 2354, SPSS 1152 and SPSS 
1501. For observations of other systems we refer the reader to 
table 1 in the following publications: Feline et al. 2004a (OU 
Vir), Feline et al. 2004b (XZ Eri and PV UMa), Littlefair 
et al. 2006a (SPSS 1702), Littlefair et al. 2007 (SPSS 1507) 
and Littlefair et al. 2008 (SPSS 0903, SPSS 1227, SPSS 
1433, SPSS 1501 and SPSS 1502). 

For reasons outlined in section 13.31 we do not use the 
SPSS 1501 data listed in Littlefair et al. (2008). Instead we 
model a single eclipse observed in 2004 (Table [1]). Not all of 
the eclipses listed in Table[T]are used for determining system 
parameters. This is because the eclipses have poor signal- 
to-noise, or lack clear bright spot features. The eclipses not 
used for determining system parameters are however still 
used to refine our orbital ephemerides (section 13. 1|) . These 
eclipses include CTCV 2354 cycle numbers 11197, 11198, 
11366, 11396, 11457, 11472 and SPSS 1501 cycles 24718 and 
24719. CTCV 2354 cycle numbers 16156, 16676 and CTCV 
1300 cycle number 12888 are analysed separately in section 
I3.4l because the shape of the eclipse has changed significantly 
in comparison to the 2007 data (see section [372]) . 

Pata reduction was carried out in a standard man- 
ner using the ultraca m pi peline reduction softw are, as de- 
scribed in lFelinel |2005h and lPhillon et al.1 (|2007f ). A nearby 
comparison star was used to correct the data for trans- 
parency variations. Observations of the standard stars G162- 
66, G27-45 and G93-48 were use d to correct the m agnitudes 
to the standard SPSS system (|Smith et all 120021 '). Pue to 
time constraints and poor weather, we were unable to ob- 
serve a standard star to flux calibrate our data for SPSS 
1152. Consequently, we have used the Sloan magnitudes of 
the comparison stars and corrected for different instrumen- 
tal response. To do this, we use measured response curves 
for filters and dichroics to create overall response curves for 
ULTRACAM. These are then combine d with curves of theoret- 
ical extinction and library spectra (|Pickles!ll998l l to obtain 
synthetic ultracam colours. The same process is then re- 
peated for the SPSS colour set, with the difference between 
the two sets being the correction applied. 



Object 


T (HJD) 


forb (d) 


CTCV 1300 
CTCV 2354 
SDSS 1152 
SDSS 1501 


2454262.599146 (8) 
2454261.883885 (5) 
2455204.601298 (6) 
2453799.710832 (3) 


0.088940717 (1) 
0.065550270 (1) 
0.067721356 (3) 
0.0568412623(2) 



3 RESULTS 

3.1 Orbital ephemerides 

The times of white dwarf mid-ingress T w , and mid-egress 
T we were determined by locating the minimum and maxi- 
mum times, respectively, of the smoothed light-curve deriva- 
tive. Mid-eclipse times, Tmid, were determined by assum- 
ing the white dwarf eclipse to be symmetric around phase 
one and taking T mi d = (T W i + T we )/2. Eclipse times were 
taken from the lite rature for C TCV 1300, CTCV 2354 
i|Tappert et al1l2004ri. SPSS 1501 dLittlefair et alj|2008h and 
SPSS 1152 (|Southworth et al.ll2010! ) and combined with our 
mid-eclipse times shown in Table Q] The errors on our data 
were adjusted to give \ 2 — 1 with respect to a linear fit. In 
each case we observe no cycle ambiguity. We do however, 
observe a significa nt, O-C offset betwe en our data and the 
times published bv lTappert et all l|2004h for CTCV 1300 and 
CTCV 2354. For CTCV 1300, the average difference is 165.9 
seconds, while for CTCV 2354 it is 148.0 seconds. We believe 
this to be due to the differing methods of calculating T m id\ 
Tappert et al. (2004) calculated T mi d by fitting a parabola 
to the overall eclipse structure, whereas we determined T m id 
from the white dwarf eclipse. We therefore subtract these av- 
erage offsets from the published literature values and take 
the resulting O-C difference as our uncertainty on that time. 
Where possible, we averaged the measured mid-eclipse times 
in the r' and g' bands to fit the ephemeris. Pue to the low 
signal-to-noise of the eclipse of CTCV 2354 in May 2010, 
measuring the times of white dwarf mid-ingress T w i and mid- 
egress T we was not possible. Consequently, the eclipse times 
and errors were measured by eye in the g' and r' bands and 
then averaged in order to estimate the cycle number. This 
time was not used to refine the ephemeris but is included for 
completeness. The ephemerides found are shown in Table [2] 



3.2 Light curve morphology and variations 

3.2.1 CTCV J 1300-3052 

Fig. [1] (top) shows the two observed eclipses of CTCV 1300 
from the 2007 data set folded on orbital phase in the g' band. 
The white dwarf ingress and egress features are clearly visi- 
ble at phases 0.965 and 1.040, respectively, as are the bright 
spot features at phases 0.960 and 1.085. These features dom- 
inate the light curve, which follows a typical dwarf nova 
eclipse shape (e.g. Littlefair et al. 2006a, 2007, 2008). The 
depth of the bright spot eclipse indicates that the bright 
spot is the dominant source of light in this system, while 
the eclipse of the accretion disc is difficult to discern by eye, 
indicating that the accretion disc contributes little light to 
this system. Closer inspection of the eclipses from each night 
reveals a noticeable difference in the shape of the bright spot 
ingress feature. This is caused by heavy pre-eclipse flickering, 
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Table 1. Journal of observations. The dead-time between exposures was 0.025 s for all observations. The relative GPS time stamping on 
each data point is accurate to 50 fis. Instr setup denotes the telescope (WHT, NTT or VLT) and instrument used for each observation, 
where UCAM and USPEC represent ULTRACAM and ULTRASPEC, respectively. Phase Cov corresponds to the phase coverage of the eclipse, 
taking the eclipse of the white dwarf as phase 1. T exp and N exp denote the exposure time, and number of exposures, respectively. 



Date 




Object 


Instr setup 


T w CHMDJ1 


Cycle 


Phase Cov 


Filters 


T„^„ (s) 

-<- exp \ a j 


exp 


Sciiii 


e (") 


2007 June 


09 


CTCV 


2354 


VLT+UCAM 


54261.383926(25) 





0.73- 


-1.08 


u 


9 


r' 


2.22 


821 


0.6- 


-1.0 


2007 June 


13 


CTCV 


2354 


VL1 + UCAM 


54265.316786(61) 


60 


0.70- 


-1.11 


u 


9 


r' 


4.92 


173 


0.8- 


-1.2 


2007 June 


15 


CTCV 


2354 


VLT+UCAM 


54267.348921(21) 


91 


0.78- 


-1.08 


u 


9 


r' 


2.22 


779 


0.6- 


-0.7 


2007 June 


15 


CTCV 


2354 


VLT+UCAM 


54267.414476(20) 


92 


0.74- 


-1.05 


a 


(J 

is 


V 


2.22 


757 


0.6- 


0.7 


2007 June 


16 


CTCV 


2354 


VLT+UCAM 


54268.397717(20) 


107 


0.72- 


-1.06 


u 


9 


r' 


2.22 


845 


0.6- 


-1.1 


2007 June 


19 


CTCV 


2354 


VLT+UCAM 


54271.413077(21) 


153 


0.82- 


-1.15 


u 


9 


r' 


2.22 


826 


0.6- 


-1.0 


2007 June 


20 


CTCV 


2354 


VLT+UCAM 


54272.396368(29) 


168 


0.50- 


-1.50 


u 


9 


r' 


2.32 


2390 


0.6- 


-1.0 


2007 June 


21 


CTCV 


2354 


VLT+UCAM 


54273.314054(5) 


182 


0.86- 


-1.20 


u 


9 


r' 


1.96 


931 


1.2 


-2.4 


2007 June 


21 


CTCV 


2354 


VLT+UCAM 


54273.379579(3) 


183 


0.77- 


-1.09 


LI 


9 




1.96 


916 


0.9- 


-1.5 


2009 June 


12 


CTCV 


2354 


NTT+USPEC 


54995.350263(6) 


11197 


0.65- 


-1.35 




9 




9.87 


182 


1.4- 


-2.6 


2009 June 


12 


CTCV 


2354 


NTT+USPEC 


54995.415961(6) 


11198 


0.35- 


-1.17 




a 




9.87 


482 


1.0- 


-2.2 


2009 June 


23 


CTCV 


2354 


NTT+USPEC 


55006.428224(2) 


11366 


0.66- 


-1.16 




9 




3.36 


817 


1.4 


-2.2 


2009 June 


25 


CTCV 
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and is clearly visible in Fig. [2] The nickering is reduced be- 
tween phases corresponding to the white dwarf ingress and 
bright spot egress, indicating the source of the flickering is 
the inner disc and/or bright spot. Due to the heavy flicker- 
ing, we decided to fit each night individually rather than fit 
to a phase-folded average, in order to provide a more robust 
estimation of our uncertainties. Our 2010 observations are 
discussed in section [3.4. II 



3.2.2 CTCV J23 54-4700 

Fig. [1] (middle) shows all of the observed eclipses from the 
2007 data set of CTCV 2354 folded on orbital phase in the g' 
band. The white dwarf ingress and egress features are clearly 
visible and along with the accretion disc dominate the shape 
of the light curve. A weak bright spot ingress feature is vis- 
ible at an orbital phase of 0.995, however the system suf- 
fers from heavy flickering, making it difficult to identify the 
bright spot egress. The shape of the average light curve in- 
dicates possible egress features at phases 1.060 and 1.080, 
but given the scatter we cannot be certain whether these 
represent genuine egress features or merely heavy flickering. 
The flickering is reduced between phases corresponding to 
the white dwarf ingress and egress, indicating the source of 
the flickering is the inner disc. Our observations from 2009 
and 2010 are discussed in section 13.4.21 



3.2.3 SDSS J1152+4049 

Fig. [T] (bottom) shows all of the observed eclipses from the 
2010 data set of SDSS 1152 folded on orbital phase in the 
g' band. The signal-to-noise ratio of our data is low in com- 
parison to other systems, but we still see a clear bright spot 
ingress feature at phase 0.975 in addition to a clear bright 
spot egress feature at phase 1.075. The white dwarf features 
are clear, and dominate the overall shape of the light curve. 
Like CTCV 1300, the eclipse of the accretion disc is difficult 
to discern by eye, which again suggests that the accretion 
disc contributes little light to this system. 

3.3 Light curve modelling 

To determine the system parameters we used a physical 
model of the binary system to calculate eclipse light curves 
for the white dwarf, bright spot, accretion disc and donor. 
Feline et al. (2004b) showed that this method gives a more 
robust determination of the system parameters in the pres- 
ence o f flickering than the derivative method of lWood et al.l 
(1986). The model itself is based on the techniques devel- 
oped by IWood et al.l (|l985l ) and iHorne et all ll 19941). and 
is an adapted version of the one used by iLittlefair et alj 
(2008). This model relies on three critical assumptions: the 
bright spot lies on the ballistic trajectory from the donor 
star, the donor fills its Roche lobe, and the white dwarf is 
accurately described by a theoretical mass-radius relation. 
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Figure 1. ULTR.ACAM g' band light curves of CTCV 1300 (2007, 
top), CTCV 2354 (2007, middle) and SDSS 1152 (2010, bottom). 



Obviously these assumptions cannot be tested directly, but 
it has been shown that the masses derived with this model 
are consistent with other methods com monly employed in 
CVs over a range of orbital periods (e.g. iFeline et al.ll2004bl; 
iTulloch et ai]|2009l; ICopperwheat et alJ l201Ch . The model 
used bv lLittlefair et al.l ([2008) had to be adapted due to the 
prominent bright spot observed in CTCV 1300 (Fig. [J). The 
'old' model fails to correctly model the bright spot ingress 
and egress features in this system satisfactorily and results 
in a poor fit. We have thus adapted the model to account 
for a more co mplex bright spot by addin g four new parame- 
ters, following ICopperwheat et al.l (|2010h . bringing the total 
number of variables to 14. These are: 



(i) The mass ratio, q = M r /M w . 

(ii) The white dwarf eclipse phase full-width at half- 
depth, A(f>. 

(hi) The outer disc radius, Rd/a, where a is the binary 
separation. 

(iv) The white dwarf limb-darkening coefficient, U w . 

(v) The white dwarf radius, R w /a. 

(vi) The bright-spot scale, S/a. The bright spot is mod- 
elled as two linear strips passing through the intersection 
of the gas stream and disc. One strip is isotropic, while the 
other beams in a given direction. Both strips occupy the 
same physical space. The intensity distribution is given by 
(X/5') y e - ' x / S ' ) , where X is the distance along the strips. 
The non-isotropic strip does not beam perpendicular to its 
surface. Instead the beaming direction is defined by two an- 
gles, Quit and 6 ya w 

(vii) The first exponent, Y, of the bright spot intensity 
distribution. 

(viii) The second exponent, Z, of the bright spot intensity 
distribution. 



(ix) The bright-spot angle, Q az , measured relative to the 
line joining the white dwarf and the secondary star. This 
allows adjustment of the phase of the orbital hump. 

(x) The tilt angle, Quit, that defines the beaming direc- 
tion of the non-isotropic strip. This angle is measured out of 
the plane of the disc, such that Quit = would beam light 
perpendicular to the plane of the disc. 

(xi) The yaw angle, Q ya w This angle also defines the 
beaming direction of the non-isotropic strip, but in the plane 
of the disc and with respect to the first strip. 

(xii) The fraction of bright spot light that is isotropic, 

fiso- 

(xiii) The disc exponent, b, describing the power law of 
the radial intensity distribution of the disc. 

(xiv) A phase offset, <f>o. 

The data are not good enough to determine the white 
dwarf limb-darkening coefficient, U w , accurately. To find an 
appropriate limb-d arkening coeffic i ent, w e follow the pro- 
cedure outlined in iLittlefair et al.l (|2007l ). whereby an es- 
timate of the white dwarf effective temperature and mass 
is obtained from a first iteration of the fitting process out- 
lined below, assuming a limb-darkening coefficient of 0.345. 
ILittlefair et al.l (|2007 ) show that typical uncertainties in U w 
are ~ 5%, which leads to uncertainties in R w /a of ~ 1%. 
These errors have negligible impact on our final system pa- 
rameters. 

As well as the parameters described above, the model 
also provides an estimate of the flux contribution from the 
white dwarf, bright spot, accretion disc and donor. The 
white dwarf temperature and distance are found by fitting 
the white dwarf fluxes from our model to the prediction s of 
white dwarf model atmospheres (|Bergeron et al.lll995l ). as 
shown in Fig. [3] We find that with the exception of CTCV 
2354, all of the systems analysed lie near, or within, the 
range of white dwarf col o urs a llowed by the atmosphere 
models of [Bergeron et all (|l995l ). although the systems do 
not always lie near the tr ack for the appr o priate mass and ra- 
dius of the white dwarf. ILittlefair et al.l (|2008l ) compare the 
temperatures derived using light curve fits to those found us- 
ing SDSS spectra and GALEX (Galaxy Evolution Explorer) 
fluxes for a small number of systems and conclude their 
white dwarf temperatures are accurate to ~1000K. The sys- 
tems examined by Littlefair et al. (2008) are all found to lie 
close to the Bergeron tracks; it is likely that systems that lie 
far from the tracks are less accurate. We no te our tempera- 
tures have larger uncertainties than those of ILittlefair et aLl 
(2008). This is because our temperatures take into account 
the uncertainty in white dwarf ma ss when compa r ing th e 
white dwarf fluxes to the models of lBergeron et akl (|l995l ). 

It is possible that our white dwarf colours are affected 
by contamination from the disc or bright spot, or an un- 
modelled light source such as a boundary layer. If our white 
dwarf colours are incorrect, then our derived white dwarf 
temperatures will be affected. Changing the white dwarf 
temperature will alter U w . Our model fitting measures R w /a 
and uses a mass-radius relationship to infer M w , which is 
then used to find the mass of the donor star. However, U m 
and R w are partially degenerate, so U w therefore affects R m 
and M w . M w is also affected by temperature changes be- 
cause the white dwarf mass-radius relationship is temper- 
ature dependent. The white dwarf temperature also affects 



6 C. D. J. Savoury et al. 



the luminosity of the system, and hence distance estimate. 
It is therefore important to quantify the effect that incor- 
rect white dwarf temperatures may have on distance esti- 
mates and our final derived system parameters. To do this, 
we altered the white dwarf temperature by 2000K and per- 
formed the fitting procedure described above on our best 
quality, white-dwarf dominated systems. For lower quality 
data, the random errors dominate over any systematic er- 
rors, and thus changes to the best quality data represent a 
worst case scenario. We find that changing the white dwarf 
temperature by 2000K changes R w /a by less than lcr. The 
white dwarf distance estimates change by 10-20pc. We there- 
fore conclude any error in white dwarf temperature that may 
occur does not affect our final system parameters by a sig- 
nificant amount. We note here that our moedelling does not 
include treatment of any boundary layer around the white 
dwarf, and assumes all of the white dwarf's surface is visi- 
ble. Either effect c ould lead to system atic uncertainty in our 
white dwarf radii l|Wood et al.lll986l ). 

A Markov Chain Monte Carlo (MCMC) analysis was 
used to adjust all parameters bar Uw- MCMC analysis is 
an ideal tool as not only does it provide a robust method for 
quantifying the uncertainties in the various system param- 
eters, it is more likely to converge on the global minimum 
\ 2 ra t her t h an a local minim um \ 2 ■ We refer the reader to 
iFordl HH), iGreeorvl (|2007l ) and references therein for ex- 
cellent overviews of MCMC chains and Bayesian statistics 
and limit ourselves to a simple overview. 

MCMC is a random walk process where at each step in 
the chain we draw a set of model parameters from a nor- 
mal, multi-variate distribution. This is governed by a co- 
variance array, which we estimate from the initial stages of 
the MCMC chain. The step is either accepted or rejected 
based on a transition probability, which is a function of the 
change in \ 2 ■ We adopt a transition probability given by the 
Metropolis-Hastings (M-H) rule, that is P = exp" Ax2/2 . 
The sizes of the steps in the MCMC chain are multiplied 
by a scale factor, tuned to keep the acceptance rate near 
0.23, which is found t o be the optimal valu e for multi-variate 
chains such as these (|Roberts et alj[l997l l. 

A typical MCMC chain included some 700,000 steps, 
split into two, 350,000 step sections. The first section is used 
to converge towards the global minimum and estimate the 
covariance matrix (known as the burn- in phase). The sec- 
ond section fine tunes the solution by sampling areas of pa- 
rameter space around the minimum. In doing so, we also 
produce a robust estimation of our uncertainties. Together, 
these steps are usually sufficient to enable the model to con- 
verge on the statistical best fit, regardless of the initial start- 
ing parameters. 

While implementing the MCMC code, we discovered a 
bug in our original code. The re-binning code used to av- 
erage several light curves together mistreated the widths of 
the bins, which in turn affected the trapezoidal integration of 
the model over these bins. The direct result was that in cases 
of heavy binning, such as systems with heavy flickering or 
where several light curves had been averaged together (e.g. 
SDSS 1502), the white dwarf radius, Rw/a, was underesti- 
mated. The exact amount depended on the level of binning 
used. This consequently resulted in an overestimate of the 
white dwarf mass. Since the mass of the donor star, M r , 
is related to the white dwarf mass M m by M r = qM w , we 



were also left with an overestimate of the donor mass. This 
problem affects all of our previously published eclipsing-CV 
papers (Feline et al. 2004a, 2004b; Littlefair et al. 2006a, 
2006b, 2007, 2008) by differing amounts. However, in most 
cases re-modelling provides new system parameters that are 
within 1 — 2a of our original results, with only two exceptions 
(see section [3^4)) . The new results are presented in Table [3] 

For each system we ran an MCMC simulation on each 
phase- folded u , g' , r' or i' light curve from an arbitrary 
starting position. Exceptions include CTCV 1300, where 
each night of observations was fit individually, and SDSS 
1152 and SDSS 1501, for which we only calculated fits in 
the g' and r' bands due to u'-band data of insufficient qual- 
ity to constrain the model. Where no v! band MCMC fit 
could be obtained, we fit and scaled the g' band model to 
the u' band light curves without \ 2 optimisation. This al- 
lows us to estimate the white dwarf flux in the it' band, 
and thus estimate the white dwarf temperature. In the case 
of SDSS 1501, we also fit a d i fferen t data set to the 2006 
WHT data of Littlefair et alJ l|2008h . We fit our model to 
the single light curve dated 2004 May 17. This 2004 data 
was not fit by Littlefair et al. (2008) as the simplex meth- 
ods used gave a seemingly good fit to the 2006 data. De- 
spite appearing to have converged to a good fit, the MCMC 
analysis revealed that the 2006 data does not constrain the 
model, most likely due to the very weak bright spot features. 
The 2004 data shows much clearer and well-defined bright- 
spot features than the 2006 data (see Fig.l of Littlefair et 
al. 2008), and so despite only having one eclipse (and thus 
lower signal-to-noise) it is favoured for the fitting process. 
In general our fits to each system are in excellent agreement 
with the light curves (see Fig. [2]), giving us confidence that 
our new models accurately describe each system. 

To obtain final system parameters we combine our 
MCMC chains with Kepler's 3 rd law, the orbital period, 
our derived white dwarf temperature, and a series of white 
d warf ma s s-radi us relationships. We favour the relationships 
of IWoodl (|l995h . because they have thicker hydrogen lay- 
ers which may be more appropriate for CVs. However, they 
do not reach high enough masses for some of our systems. 
Above M w = l.OMp ) , we adopt the mass-radius relation- 
ships of iPanei et all <|2000t ) . In turn, these models do not 
extend beyond M w = 1.2 Ma; above this mass we use the 
lHamada Sz Salpeterl f|l96lh relationship. No attempt is made 
to remove discontinuities from the resulting mass-radius re- 
lationship. 

We calculate the mass ratio q, white dwarf mass 
M w /Mq, white dwarf radius R w /Rq, donor mass M t /Mq, 
donor radius R t /Rq, inclination i, binary separation a/R® 
and radial velocities of the white dwarf and donor star (K m 
and K r , respectively) for each step of the MCMC chain. 
Since each step of the MCMC has already been accepted or 
rejected based upon the Metropolis-Hastings rule, the distri- 
bution function for each parameter gives an estimate of the 
probability density function (PDF) of that parameter, given 
the constraints of our eclipse data. We can then combine the 
PDFs obtained in each band fit into the total PDF for each 
system, as shown in Fig. [4] We note that most systems have 
system parameters with a Gaussian distribution with very 
little asymmetry. Our adopted value for a given parameter 
is taken from the peak of the PDF. Upper and lower error 
bounds are derived from the 67% confidence levels. For sim- 
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plicity, since the distributions are mostly symmetrical, we 
take an average of the upper and lower error bounds. The 
final adopted system parameters are shown in Table O al- 
though Fig. [7| and [8] show the true 67% confidence levels for 
the white dwarf mass and donor mass respectively, for each 
system. 

3.4 Notes on individual systems 

3.4.1 CTCV 1300 

As noted in section [331 the two eclipses of CTCV 1300 were 
modelled individually due to the heavy pre-eclipse flicker- 
ing present in each light curve. The u'g'r' eclipses from the 
night of 2007 June 10 gave consistent results, as did the 
u'g'i' eclipses from the night of 2007 June 13. However, re- 
sults from the individual nights were not consistent with 
each other, and we are presented with two distinct solutions 
which are the result of heavy pre-eclipse flickering altering 
the shape of the bright spot ingress feature. This in turn 
gives two very different values for the mass ratio, q. Encour- 
agingly we find our white dwarf masses and radii are consis- 
tent between nights. To derive the final system parameters 
we use the PDFs as outlined in section [3731 and then take an 
average of the solution from each night. The error is taken 
as the standard deviation between the two values. 

In Fig. [5] we see a (/-band eclipse from June 2010 to- 
gether with a modified version of the model obtained from 
the 2007 dataset. This new model is found by starting from 
an average of the two 2007 fits and using a downhill sim- 
plex method to vary all parameters bar q, A<j>, R m /a and 
U w . These parameters should not change with time, and so 
the simplex fit will confirm if our bright spot positions and 
white dwarf dwarf radius are correct, and thus if our system 
parameters are reliable. We use a simplex method for two 
reasons; firstly since we are performing a consistency check 
and are not extracting system parameters from the fitting 
process we do not require a full MCMC analysis. Secondly, 
the bright spot flux appears to have reduced significantly, 
and the strength of the bright spot ingress feature means 
that we cannot constrain a full fit using the MCMC model 
used previously. The model confirms the bright spot flux 
has decreased considerably, although an orbital hump is still 
visible. The white dwarf flux remains almost unchanged, al- 
though the disc appears brighter. The fit to the data is good, 
indicating that the models derived from the 2007 data (and 
used to derive our system parameters) are reliable. 

3.4.2 CTCV 2354 

In section [3. 2. 2 1 we noted that the shape of the light curve in 
Fig. [l] indicated possible bright spot egress features around 
phases 1.060 and 1.080. Fig. [2] shows that our model has 
fit the bright spot egress feature at phase 1.080. Given the 
strength and shape of the bright spot features and general 
scatter present in the light curve we cannot be certain if the 
bright spot positions have been correctly identified by our 
model, and thus there is some element of doubt as to the 
value obtained for our mass ratio and thus donor mass. It is 
at this point we draw the readers attention to our 2010 data, 
shown in Fig. [6] The eclipse dated 2010 May 3 (centre panel) 
shows clear bright spot ingress and egress features, with a 
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Figure 5. Our June 2010 g' band eclipse of CTCV 1300, together 
with a modified model found using the 2007 data. Starting from 
an average of the two 2007 models, we ran a downhill simplex fit 
varying all parameters bar q, A0, R w /a and U w . 

clear orbital hump visible from phases 0.70-0.95. The system 
is much brighter in this state than the 2007 data previously 
modelled, in part due to a dramatic increase in bright spot 
flux. As with CTCV 1300 we carried out a downhill simplex 
method to vary all parameters bar q, A<j), R w /a and U w . 
Our fit to the eclipse of May 3 is especially pleasing, as its 
excellent agreement with the light curve confirms that our 
2007 model correctly identified the bright spot egress feature 
and thus the mass ratio obtained is reliable. 

Fig. [5] also shows a single eclipse observed just one 
month later (June 2010, right panel), and six eclipses av- 
eraged together from June 2009 (left panel). Both of these 
datasets are fit with a downhill simplex model as above. The 
June 2009 and June 2010 datasets are in stark contrast to 
the May 2010 data, with the bright spot features appear- 
ing extremely faint (2009), or seemingly non-existent (June 
2010). The disc flux in the June 2010 data appears to have 
increased significantly, giving rise to a distinct "u" shape. 
The rapid change in bright spot and disc light curves over 
such short (1 month - 1 year) time scales suggests that the 
disc is highly unstable. 

3.4.3 DV UMa 

Our donor mass derived for DV UMa has increased by 
6.1ct (AM=O.O39M ) from that published by Feline et al. 
(2004b). Close inspection of the original fit reveals that the 
bright spot features are fit poorly by the old bright spot 
model. This arises because the old bright spot model could 
not describe the complex bright spot profile present and an 
innacurate value of the mass ratio is found as a result. Our 
new bright spot model is much better in this respect, and 
is able to take into account a wider variety of geometric ef- 
fects and orientations. Given that our white dwarf radius is 
consistent with that of Feline et al. (2004b), this seems the 
most likely cause of such a large change. It is worth noting 
that our new donor masses for both DV UMa and XZ Eri, 
are both consistent with the masses obtained by Feline et al. 
(2004b) using the derivative method, which, unlike our pa- 
rameterised model, does not make any attempt to recreate 
the bright spot eclipse profile (e.g. Wood et al. 1986; Home 
et al. 1994; Feline et al. 2004a; Feline et al. 2004b). 
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Table 3 (system parameters) goes here (landscape). See separate file (tablc3params) . 

Table 3. 
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Figure 2. The phased-foldcd u'g'r 1 or u'g'i' light curves of the CVs listed in Tabic [3] fitted using the model outlined in section [3.31 
The data (black) are shown with the fit (red) overlaid and the residuals plotted below (black). Below are the separate light curves of the 
white dwarf (dark blue), bright spot (light blue), accretion disc (green) and the secondary star (purple). Data points omitted from the 
fit are shown in red. \ 2 values for each fit, together with the number of degrees of freedom (DOF) are also shown. 



3.44 SDSS 1502 

Our new fits to SDSS 1502 decrea se the donor mass b y 
2.9a (AM r = O.O12M ) from that of lLittlefair et all |200St ). 
Our mass ratio and i nclination are consistent with those of 
iLittlefair et all l^OOct ). however our white dwarf radius, R w , 
has increased by 13 percent (3.4<r). We believe the primary 



reason for this change was that the original fit was heavily 
binned, and thus more susceptible to the bug outlined in 
section [3.31 
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Figure 2 - continued 



3.4.5 SDSS 1501 

The most important change of all of our re-modelled sys- 
tems is for that of SDSS 1501. Whi lst our donor mass ha s 
only increased by 1.9a from that of iLittlefair et af. (2008), 
we note that our uncertainties are large (aM r = O.OIOM©) 
and the mass difference is large enough to take this system 
from being a post-period-bounce system, to a pre-period- 
bounce system. Although our errors do not formally rule 



out the possibility of SDSS 1501 being post-period-bounce, 
the donor's position with respect to the evolutionary tracks 
shown in Fig.[8]strongly favours that of a pre-period-bounce 
system. Such a change arises from a differenc e in bright spot 
positi ons between our model and that of ILittlefair et al.l 
l|2008h . whic h in turn affect the m ass ratio obtained. The 
data used bv lLittlefair et ah I (|2008l ) shows a very weak bright 
spot ingress feature. With the improvements made to the 
modelling process resulting from the introduction of MCMC, 



Short Period CVs below the Period Gap 11 







2 = 350, DOF, 288 " 
(WlLLLLi 










SDSS 0903 U' 



X = 640, DOF: 288 ■ ° 




SDSS 0903 



1 | 



0.95 


1 1.05 




— i— i— | — i— i— i— i — |— i— i— i — i | i i — i— i— 

X z =139, DOF: 201 


SDSS 1227 

_J 1 1 1 1 1 1 1 L 


xC 


0.9 0.85 


1 1.05 1.1 



>> d 
E 





s 


X : = 235,DOF: 158 








... .n 1 .... 1 l.MkJ 


SDSS 1433 


u' 


0.9 0.95 1 1.05 1.1 



0.S5 1 


1.05 




X ~ 589, DOF: 406 












: 


'. SDSS 1227 




j , i i i i i i i i i i_ 


0.9 0.95 1 


1.05 1.1 




— r— I— I— 1 — I — I— 



X - 556, DOF: 279 



i i i i i i 



0.9 0.95 1 1.05 1.1 




5 


-i i i i i i i i i i i i 

X 2 = 652, DOF: 308 " 

3n il . DlISi 1 






[ SDSS 0903 f 

i i i i i i i i . — i — . — i — i — i — i — . — . — . — . — 



0.95 


1 1.05 




X 2 = 478, DOF: 398 \ 






— ! 




SDSS 1227 


T* 


0.9 0.95 


1 1.05 1.1 







X 2 =" 443. DOF: 276 

[ 












SDSS 1433 

.... i .... i . . 


r* 

. . 1 . 




O.S 0.95 


1 


1.05 1.1 





X 


= 2948, DOF: 311 






"»«IWv^* 

SDSS 1502 


r' 


0.9 0.95 1 


1.05 1.1 




( 1 7/05/04) 



0.9 0.95 1 1.05 1.1 

Orbital Phase 



0.9 0.95 1 1.05 1.1 

Orbital Phase 



Figure 2 - continued 



it is clear that the 2006 data used bv iLittlefair et all (|2008l ) 
does not constrain the mass ratio, q, tightly enough. In con- 
trast, the 2004 data shows much clearer bright spot features, 
and we therefore favour this over the 2006 data as discussed 
in section [3731 



4 DISCUSSION 

4.1 White dwarf masses 

Population studies bv lWillems et al l (|2005l ) predict that be- 
tween 40 and 80 percent of CVs are born with He-core 
white dwarfs (M w < O.5OM0) and therefore He-core white 
dwarfs (He-WDs) are expected to be common amongst CV 
primaries. It is surprising then that out of our sample of 
14 systems, we observe no He-WDs. Of all of our objects, 
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Figure 3. The white dwarf colours derived from our model fitted together with the white dwarf models 


of iBereeron et al.l (Il995ll 


top to bottom, each 


curve represents log g = 9.0, 8 


5, 8.0, 7.5 and 7.0 respectively. The measured white dwarf colours 


are shown i 



red, and are used to derive the white dwarf temperature, which in turn is used to correct the white dwarf mass-radius relationships used 
later to obtain the final system parameters. 



SDSS 1152 is found to have to the lowest white dwarf 
mass with M w = .560 ± 0.028. The mass distribution of 
iKepler et al. | (|2007h for SDSS white dwarfs suggests He- 
WDs have a typical mass of ~O.38M0. The most massive 
He-WDs are thought to form from single RGB stars, which 
due to extreme mas s loss are able to avoid the Helium flash; 
iD'Cruz et al.l |l996) consider models with a range of mass 



loss rates on the RGB and manage to produce He-WDs with 
masses up to ~O.48M0. It is likely that this represents an 
upper limit to the mass of He-WDs and hence SDSS 1152 
is too massive to be a viable candidate for a He-core white 
dwarf. 

We find that our white dwarf masses are not only too 
massive to be He-WDs, but are also well in excess of the aver- 
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ag e mass for singl e DA white dwarfs. Using the same method 
as iKniggel (|2006l ). we calculated the average white dwarf 
mass of our entire sample to be M w = 0.81 ± 0.04Mp), with 
an int rinsic scatter of O.13M0. In comparison. iLiebert et all 
|2005l ) find the mean mass of DA white dwarfs to be 
M w ~0.603 M Q , while lKepler et alJ l|2007l ) find a mean mass 
of M w = 0.593 ±O.O16M . 

Ou r stu d y thus supports p revious findings (e.g . 

1976; 



Uu r stu d y tnus supports p revious nndmgs 
IWarnerl Il973l . Il976l ; iRitterl 1 19761 . 1 19851 ; iRobinsonl 



ISmith fc Dhillonlll998l ; IKniggel [2006l 'l that white dwarfs in 
CVs are o n average much highe r in mass than single field 
stars. Like lLittlefair et alJ l)2008h . we compare our masses to 
the average mass of M w ~ 0.7 3 ± 0.05Mp) fo r white dwarfs 
in CVs below the period gap (jKniggc 2006) and find that 
our white dwarf masses are generally much higher. This is 
especially so for systems P or t ^ 95 mins, where we find a 
mean white dwarf mass of M m — 0.83 ± O.O2M0, with an 
intrinsic scatter of O.O7M0. Some our white dwarf masses 
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Figure 6. Our g' band observations of CTCV 2354 from 2009 (left), 2010 May 3 (centre) and 2010 June 6 (right), fit with a modified 
version of the model obtained from the 2007 data as described in section 13. 1.21 



are revised, and have accordingly m oved downwards in mass 
compared to Littlefair et ail {2008), but as Fig. Qshows, for 
Porb ^ 95 mins, 8 out of 9 sy s tems are more massive than 
the average found by IKniggd {20061 ), with SDSS 1502 the 
only exception. On the s ame plot, we a lso plot the disper- 
sion of masses as found bv lKniggel {20061) and the mean mass 
of sing le SDS S white dwarfs found by iKepler et~aT1 {20o3). 
Note that the IKniggd {2006) sample contains three systems 
also included in our study: OU Vir, DV UMa and XZ Eri, 
using the old mass determinations of Feline et al. (2004a, 
2004b). We see most of our masses are within the dispersion 
found, indicating that no individual white dwarf mass is un- 
usual. However 8 out of 9 systems above average does seem 



anomalously high, considering that if we model the white 
dwarf masses as a Gaussian distribution, the probability of 
such an occurrence is less than 2 percent (independent of 
the actual mean or variance). Suc h difference between our 
sample and that of Knigge (2006) is concerning, and it is 
therefore desirable to consider the selection effects, consid- 
ering th e majority of our short period system s are all SDSS 
objects {Szkodv et al.ll2004l . l2005l . l2006l , l2007r i. 

The majority of SDSS CVs found are generally rejected 
quasar candidates with a limiting magnitude of g' = 19— 20, 
and are initially selected for fol l ow up on the basis of u'-g' 
colour cuts ( see iGansicke et al.l (l200Sh for a more in depth 
description) . iLittlefair et al. ( 20081 ) show that systems with 
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M w ^ O.5OM0 are blue enough to pass the SDSS colour cuts 
and conclude that selection effects such as these are unlikely 
to explain the high mass bias of our white dwarf sample. 
However, the majority of the systems we have studied are 
close (g'~ 17.5 - 19.5) to the g' = 19 limit of the SDSS sur- 
vey. This raises the possibility that the SDSS sample only 
finds t he brightest of the short period CVs. lRitter fc Burkertl 
1 19861 ) have shown that CVs with high mass white dwarfs 
are brighter than their low mass counterparts. This suggests 
there maybe some bias towards high mass dwarfs. However, 
this conclusion is not appro priate to the systems studied 
here: lRitter fe Burkertl (|l986T ) only consider the effects of ac- 
cretion luminosity whereas in most of our systems the white 
dwar f considerably outshin es the accretion disc. 

IZorotovic et alj (|201ll ) consider selection effects in white 
dwarf dominated SDSS systems with a variety of differ- 
ent mass transfer rates and conclude there is actually bias 
against high mass white dwarfs. They find a O.9OM0 white 
dwarf is approximately 0.15 magnitudes fainter than a 
O.75M0 white dwarf, which corresponds to a decreased de- 
tection efficiency of ~ 20%. This suggests that our finding of 
high mass white dwarfs at short orbital periods is not due to 
selection effects, and is in fact a true representation of the 
intrinsic mass distribution of CVs. However, this analysis 
only considers white dwarf luminosity; in some of our sys- 
tems the bright spot features are prominent, and contribute 
significantly to the overall flux of the system. It remains 
possible that the finding of very high white dwarf masses 
for our short period CVs is due to selection effects. How- 
ever, a full and thorough quantification of any bias would 
require detailed calculations of the luminosity of white-dwarf 
dominated systems (including bright spot emission) plus an 
investigation of the selection effects in the SDSS sample. 
Such an analysis is beyond the scope of this paper. 

Our results have important consequences for the mod- 
elling of nova outbursts and their impact on the long-term 
evolution on CVs. Typical calculations show that the mass 
of the white dwarf decre ases by between 1 and 5 percent 
per 1 000 nova cycles fe.g. lYaron et al] |2005; Epclstai n et al.l 
120071 ). The dominance of high- mass white dwarfs in our sam- 
ple of short period systems suggests that any white dwarf 
erosion due to nova explosions must be minimal, or that 
not all of the accreted matter is ejected during nova ig- 
nition, resulting in the white dwarf mass increasing over 
time. This could, in principle enable the white dwarfs in 
cataclysmic variables to grow in mass until they reach the 
Chandrasekhar limit. 



4.2 Period bounce 

Population synthesis models fo r cataclysmic variables (e.g. 
lKolblll993l ; IWillems et ai1l2005l ) all predict that large num- 
bers of the CV population (~15 - 70 percent) have evolved 
past the period minimum. This has always been in stark 
contrast to observations, possib ly in part due to selection 

effe cts (e.g. ILittlefair et alJ|2003h . 

iLittlefair et al.l (|2006bl , 120071 . 120081 ) identified four sys- 
tems (SDSS 1035, SDSS 1507, SDSS 1433, SDSS 1501) with 
donors below the sub-stellar limit, three of which are likely 
to be post-period-bounce CVs (SDSS 1035, SDSS 1501 and 
SDSS 1433). Our subsequent re-analysis gives three systems 
(SDSS 1035, SDSS 1507, SDSS 1433) with donors below 



the sub-stellar limit, two of which (for reasons outlined be- 
low) we believe are post-period-bounce CVs (SDSS 1035 and 
SDSS 1433). SDSS 1501, which no longer features as a post- 
period-bounce system, is discussed in greater detail in sec- 
tion [3T431 

ISirotkin fc Kiml (|2010l ) claim that SDSS 1433 cannot 
be considered a post-period-bounce object since the mass 
transfer rates and donor star temperatures implied are too 
high. The mass transfer rate i s found using an estimate o f 
the white dwarf temperature (|Townslev fc G ansickc 2009), 
while the donor star temperature is inferred using a semi- 
empirical relationship that is also dependent on the white 
dwarf temperat u re. Th e white dwarf temp erature used by 
ISirotkin fc Kiml (|2010l ) is that derived by ILittlefair et al.l 
(2008) from model fitting. We believe, at least in this case, 
that using M and T2 is an unreliable test of the evolutionary 
status of CVs donors, since accurate determinations of the 
white dwarf temperature are difficult to obtain. Of all the 
system parameters we have derived, the white dwarf tem- 
peratures are the least well constrained, and this does not 
take into account systematic errors. Since the white dwarf 
temperature is found using the flux from just three colours, 
and our model does not include all possible sources of lumi- 
nosity (e.g. a boundary layer), there is a good chance our 
white dwarf temperatures are affected by systematic errors 
at some level, as discussed in section [331 Instead, we focus 
on the donor star mass, M r . 

If the angular momentum loss rate is similar for sys- 
tems with identical system parameters, we expect all CVs 
to follow very similar evolutionary tracks with a single lo- 
cus in the mass-period relationship (and by analogy, mass- 
radius relationship) for CV donors, as shown in Fig. [8] 
Th e empir ic al do nor star mass-radius relationship derived 
by iKniggel (|2006l ) shows that a single evolutionary track 
does very well at describing the observed M2 — Port relation- 
ship, although the shape of that relationship is poorly con- 
strained at low masses. A single evolutionary path also ex- 
plains the presence of the "period spike" , a long sought after 
fea ture in the orbi t al per iod distribution recently identified 
by iGansicke et aTl |2009). We therefore expect there to be 
a unique donor mass corresponding to the minimum orbital 
period, below which an object becomes a period-bouncer. 
The exact mass at which this occurs is very uncertain, 
and does not ne cessarily correspond to the s ub-stellar limi t 
(|Pattersonll2009l ) . From the empirical work of lKni ggc ( 2006), 
the best estimate for M bounce is M r = 0.063 ± 0.009M©. 
Three of our systems (SDSS 1035, SDSS 1433 and SDSS 
1507) fall well below this value, although SDSS 1507 is an 
u nusual system, and is d iscussed in the following section. As 
111 ILittlefair etall (|2008l ). we do not include it in our sample 
of post-period minimum CVs. We therefore have two strong 
candidates for post-period minimum CVs (SDSS 1035 and 
SDSS 1433) from our total sample of 14 CVs (nine of which 
are SDSS systems). From this, we estimate that 14 ± 7 per- 
cent of all CVs below the period gap, and 22 ± 11 percent 
of all short period CVs (P or b ^ 95 mins) have evolved past 
the period minimum. These findings are consistent, albeit 
to a crude approximation given our small sample of objects, 
with current population synthesis models. Since all of our 
short period syste ms are SDSS CV s , we c annot rule out se- 
lection effects, but iGimsicke et~aH (|2009h have shown that 
the number of period minimum CVs found within the SDSS 
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is broadly consistent with other surveys, allowing for nor- 
malisation of survey volumes. 

4.3 SDSS 1507 

The orbital period of SDSS 1507 is far below the well- 
defined period minimum and thus the nature of this system 
is of great interest to theorists and observers. It is possible 
that this system represents the true orb ital period minimum 
as predicted by iKolb fc Barafhj (|l999l ). However, if this is 
indeed the case, we would expect a large number of sys- 
tems between orbital periods of 67 m inutes, and 83 minute s 
where the period spike is observed l|Gansicke et al.| [2009). 
These systems are not observed, and hence it is like l y that 
some other mechanism is responsible. iLittlefair et~all (I2007T I 
speculate that this system was either form ed directly from 
a whi te dwarf /brown dwarf binary, while iPatterson et al.l 
(2008) argue that the system could be a member of the halo. 
Both derive system parameters, and both obtain distance 
estimates to the system. 

Our derived system par amet ers are consistent wit h 
those of Litt lefair et ail l|2007l ) and IPatterson et all l|2008h . 
within uncertai nties. Our distance es timate is in excellent 
agreement with ILittlefair et al.l (|2007l ). which is not supris- 
ing since we both calculate the distance using the same 
methods and dataset. However, our new d istance estimate 
still p l aces the system nearer than that of IPatterson et all 
l|2008h . IPatterson et al.l (|2008l ) obtain a lower limit to the 
distance using parallax. The parallax value implies a dis- 
tance d > 175 pc, which taken alone, is consistent with 



our estimate of d = 168 ± 12 pc. Patterson combines his 
parallax with a range of other observational constraints us- 
ing Bayesian methods to yield a final distance estimate of 
d — 230 ± 40 pc. If our distance of d = 168 ± 12 pc is nearer 
the true distance, then combining with Patterson's proper 
motion measurement of 0.16"/yr yields a transverse velocity 
of d = 128 ± 9 kms -1 . This lower transverse velocity is still 
very much an outlier in the distri bution of 354 CVs shown 
in Fig.l of IPatterson et all (120081 ). Therefore, regardless of 
which distance is correct, the proper motion of SDSS 1507 
still supports halo membership. 

4.4 Exploring the standard model of CV evolution 

Fig. [8] shows the evolutionary models of IKolb tk Baraffd 
(1999) calculated with enhanced mass-transfer rates. Also 
shown is a model with 50 percent star spot coverage on the 
surface of the donor. P ositions of the period minimum, and 
period gap as found by iKnigge] (|2006l ) are also shown. Mass 
determinations for all systems presented here are included. 
We see that the standard theoretical mode ls are a poor fit 
to the data. For a given mass, the models of lKolb fc Baraffel 
l|l999l ) significantly underestimate the orbital period, and 
thus the donor radii. 

Models with enhanced mass transfer rates and star spot 
coverage do rather better at reproducing the observed donor 
masses, although the general scatter of short period systems 
makes choosing be tween these difficult. T his is in line with 
the conclusions of ILittlefair et all (|2008l ). The models be- 
gin to diverge significantly at orbital periods greater than 
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100 minutes. Unfortunately, in this regime there are few sys- 
tems with precisely known donor masses. Clearly, we require 
more mass determinations for systems with orbital periods 
between 100 and 130 minutes. 



5 CONCLUSIONS 

We present high-speed, three-colour photometry of a sample 
of 14 eclipsing CVs. Of these CVs, nine are short period 
(Porb ^ 95 minutes), and one is within the period gap. For 
each of the 14 objects we determine the system parameters 
by fitting a physical model of the binary to the observed 
light curve by \ 2 minimisation. We find that two of our nine 
short period systems appear to have evolved past the period 
minimum, and thus supports various assertions that between 
15 and 70 per cent of the CV population has evolved past the 
orbital period minimum. The donor star masses and radii are 
not consistent with model predictions, with the majority 
of donor stars being ~10 per cent larger than predicted. 
Our derived masses and radii show that this can explained 
by either enhancing themass transfer rate or modifying the 
stellar physics of the donor star to take into account star spot 
coverage. Unfortunately, we still lack enough precise donor 
masses between orbital periods of 100 and 130 minutes to 
choose between these alternatives. 

Finally, we find that the white dwarfs in our sample 
show a strong tendency towards high masses. The high mass 
dominance implies that the white dwarfs in CVs are not 
significantly eroded by nova outbursts, and may actually 
increase over several nova cycles. We find no evidence for 
He-core white dwarfs within our sample, despite predictions 
that between 40 and 80 percent of short period CVs should 
contain He-core white dwarfs. 
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Figure 8. The M2 — P or b relationship for our dataset. Mass determinations for all systems using ULTR ACAM data are included . Our three 
new systems, in addition to other objects of particular interest are labelled. The evolutionary models of lKolb fe Bar affe (1999) calculated 
with different mass- transfer rates are shown with red (dashed) and blue (dot-dashed) lines. A model with 50 percent star spot coverage 
on the su rface of the do nor is shown with an orange (dotted) line. The solid (black) line sh ows the empiri cal mass-radius relationship as 
found by Knigec (200jJ). The position of the period minimum and period gap, as found bvlKniggc (200^), are also shown. 
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